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Abstract. The motion of an extended body up to the quadrupolar structure 
is studied in the Schwarzschild background foUowring Dixon's model and within 
certain restrictions (constant frame components for the spin and the quadrupole 
tensor, center of mass moving along a circular orbit, etc.). We find a number of 
interesting situations in which deviations from the geodesic motion, due to the 
internal structure of the particle, can originate measurable effects. However, the 
standard clock-effect for a pair co/counter-rotating bodies spinning up/down is 
not modified by the quadrupolar structure of the particle. 



PACS number: 04.20.Cv 

1. Introduction 

The equations of motion for an extended body in a given gravitational background were 
deduced by Dixon [U [21 [3l HI [5] (hereafter the "relativistic extended body model," or 
simply "Dixon's model") in multipole approximation to any order. In the quadrupole 
approximation they read 

D P'^ 1 1 

= - -R^'^c^aW'S"^ - -r'^^^Rc.B-yS'^ = i^(''pi")f , (1.1) 

CLTu 2 6 

n '^fJ-^ A 

= 2P[^C/'^1 + -J°''^^^^'R''\c.t3 , (1.2) 

dT(7 3 

where P'' = m[/^ (with Up ■ Up — —1) is the total four-momentum of the particle, 
and 5''^'' is a (antisymmetric) spin tensor; U is the timelike unit tangent vector of 
the "center of mass line" Cjj used to make the multipole reduction, parametrized by 
the proper time tu. The tensor J'^P'f^ is the quadrupole moment of the stress-energy 
tensor of the body, and has the same algebraic symmetries as the Riemann tensor. 
Using standard spacetime splitting techniques it can be reduced to the following form 

where Q"^ = Q'^°^^^ represents the quadrupole moment of the mass distribution 
as measured by an observer with 4- velocity u. Similarly (with the 
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additional property ttI^'^tI = q) and n"'37<5 ^ ii^ai3]M ^^^.g essentially the body's 
momentum and stress quadrupole. Moreover the various fields Q"'^, 7r"'37 and W^f^^^ 
are all spatial (i.e. give zero after any contraction by u) . The number of independent 
components of J°'P'^^ is 20: 6 in Q"'^, 6 in Y]°'P^^ and 8 in tt"'''''. When the observer 
u = Up, i.e. in the frame associated with the momentum of the particle, the tensors 
Q"^ , TT°'^'^ and n"''''''^ have an intrinsic meaning. 

There are no evolution equations for the quadrupole as well as higher multipoles 
as a consequence of the Dixon's construction, so their evolution is completely free, 
depending only on the considered body. Therefore the system of equations is not self- 
consistent, and one must assume that all unspecified quantities are known as intrinsic 
properties of the matter under consideration. 

In order the model to be mathematically correct the following additional condition 
should be imposed to the spin tensor: 

Sf'^Up, = 0. (1.4) 

Such supplementary conditions (or Tulczyjew-Dixon conditions [6l[T]) are necessary 
to ensure the correct definition of the various multipolar terms. 

Dixon's model for structured particles originated to complete and give a rigorous 
mathematical support to the previously introduced Mathisson-Papapetrou model 
[3 [HI [9l llOj , i.e. a multipole approximation to any order which includes evolution 
equations along the "center line" for all the various structural quantities. The models 
are then different and a comparison between the two is possible at the dipolar order 
but not once the involved order is the quadrupole. 

In this paper we limit our considerations to Dixon's model under the further 
simplifying assumption [11] [12] that the only contribution to the complete quadrupole 
moment J^^T'S stems from the mass quadrupole moment Q""^ , so that tt"^'< = = 
n"^T* and 

^ -3C/^°Q^l[TC/p^l , Q^'^Upp = ; (1.5) 

Let us introduce the spin vector by spatial (with respect to Up) duality 

= ^vJ'''U^S^s , (1.6) 

where rjap^s — ^/—g^ap-^s is the unit volume 4-form and eap-yS (eoi23 = 1) is the 
Levi-Civita alternating symbol, as well as the scalar invariant 

= Is^uS''" ■ (1.7) 

In general s is not constant along the trajectory of a spinning particle. 

The assumption that the particle under consideration is a test particle means that 
its mass, its spin as well as its quadrupole moments must all be small enough not to 
contribute significantly to the background metric. Otherwise, backreaction must be 
taken into account. 



2. Motion of extended bodies in Schwarzschild spacetime 

Let us consider the case of the Schwarzschild spacetime, with the metric written in 
standard Schwarzschild coordinates. 



2M^ ^ 



ds" = - ( 1 - ^ 1 dt^ + ( 1 - — ) dr" + r^ide' + sin' ), (2.1) 
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and introduce an orthonormal frame adapted to the static observers 

ei - (1 - 2M/r)-^/^dt, e,-. = (1 - 2M/ry/^dr, = ^de, e^ = ^7^5^, (2.2) 
with dual 

J = (1 - 2M/ry^^dt, = (1 - 2M/r)-^/^dr, J = rd0, lu^ ^rsin0d(t> . (2.3) 

Let us assume that U is tangent to a (tiniehke) spatially circular orbit, hereafter 
denoted as the L/-orbit, with 

U^m+Cd^]=^[ei + iye^l ^ = {1 ~ (2.4) 

where C is the angular velocity with respect to infinity and F is a normalization factor 

r = (-ff«-C'ff00)"'^' (2.5) 

which assures that U ■ U = —1; here dot means scalar product with respect to the 
metric (j2.ip . The angular velocity ( is related to the local proper linear velocity ly 
measured in the frame (|2.2|) by 



C^J-—!^- (2.6) 



Here ( and therefore also v are assumed to be constant along the /7-orbit. We limit 
our analysis to the equatorial plane {6 = 7r/2) of the Schwarzschild solution; as a 
convention, the physical (orthonormal) component along —dg, perpendicular to the 
equatorial plane will be referred to as along the positive z-axis and will be indicated 
by z, when necessary. 

Among the circular orbits particular attention is devoted to the co- rotating (C+) 
and counter-rotating {C-) timelikc circular geodesies, having respectively C.± = ±C.k = 
±(M/r3)i/2^ so that 



M 



2M 



1/2 



2M 



3M 



1/2 



(2.7) 



with the timelike condition < 1 satisfied if r > 3Af . It is convenient to introduce 
also the Lie relative curvature pSf [14] of each orbit 

1 / 2My/^ Ck 



fc(iic) = -df In = — 1 = . (2.8) 

Let P — mUp such that 

Up = Ip [ei + vpe^] , 7p = (1 - ^ir^'^ , (2.9) 

i.e. let us assume that Up also is tangent to a circular orbit and set up an orthonormal 
frame adapted to Up given by 

eo = C/p , ei = , 62 = 7p [t-pC^ + e^] , 63 ; (2.10) 

hereafter all frame components of the various fields are meant to be referred to such a 
frame. Note that the assumption of having both U and Up aligned with a circular orbit 
certainly represents a restriction to our analysis, but leads to great simplifications 
to Dixon's equations. On the other hand, dealing with the model in its complete 
generality, i.e. with Up having no relation a priori with U, is a hard task to pursue 
analytically due to its mathematical complexity, so that a numerical analysis would 
be needed. 
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The spin vector is orthogonal to Up, so that 

S ^ S^ei+ 3^62 + S^es ■ (2.11) 
Furthermore we have 

Ooo = Qoi = Q02 = 003 = . (2.12) 

We also assume that 5*^ =0 = 5'^, and that the remaining component = S'^ = s 



(mimic of a particle moving on a circular orbit and spinning around the z-axis) as well 
as the surviving components of the mass quadrupole moment are all constant along 
the path. The latter assumption corresponds to the definition of "quasirigid motion" 
(or "quasirigid bodies") due to Ehlers and Rudolph [T^]. Clearly in a more realistic 
situation the latter hypothesis should be released. 

Consider first the evolution equations (|1.2p for the spin tensor. They imply that 

Qi2 = Qi3 = Q23 = , (2.13) 

and introducing the following "structure functions" of the extended body 

On = O33 + /, Q22 = Q33 + f, (2.14) 

they also give 

= {i^i^p - 4)s + m^{iy - Vp) + 3ukCk^ f ■ (2.15) 
Consider then the equations of motion (jl.ip . They imply that 

= i^kCk{2i^p + - f ^[-/' + 7p(2 + lyf,)!] + m{vvp - ul) . (2.16) 
2 77p 

After eliminating s through Eq. (j2.15p . the previous equation becomes 
{i^i^p - vl){~f' + 27^/) + -il[2vvpvl + vlivvp + 4) + 2vlvl]f 

2 m^Jp / 22ir,22o 2 22|4\ /r, t7\ 

~ q— + ^'^p'^K - ivVpVK -V ^k + ^k) ^ (2-17) 

where 2'^'* and 4"*^ order polynomial expressions in the velocities have been collected. 
Solving the last two equations for v and Up in terms of s and /, /' completely 
determines the motion. 

The quadrupole moment tensor of a mass distribution of density is defined 
classically by 

0(ciass) = M (ax'^x" - T^5'^'')d^x a = 1, 2, 3 (2.18) 

with = 5abX°'x^ and it is tracefree. It is therefore natural to assume the same 
property holding for the relativistic quadrupole moment tensor, obtaining 

= On + O22 + O33 = O33 + / + O33 + /' + O33 = 3033 + / + /', (2.19) 

so that the components Oafc fii this case are completely determined by the two 
"structure functions" / and /' 

Qll^\f-\f'. Q22^-\f+\f\ 033 = -^(/ + /') • (2.20) 

If the body is axially symmetric about the z-axis, then f — f and the frame 
components of Q reduce to 

Qab = diag [//3, //3, -2//3] . (2.21) 
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Without entering the problems of a relativistic definition of the quadrupole moment 
tensor Qab generahzing Eq. (|2.18p in terms of integrals over the volume of the body 
itself (the correct procedure is outlined in Dixon's works [2 [H [31 IH [S] ) , we note that 
Qab (only defined all along the world line with tangent vector U) can be interpreted as 
the the mass quadrupole moment of the extended body as measured by the observer 
Up. 

It is quite natural to introduce the following rescaled dimensionless angular and 
quadrupolar momentum quantities 

<r^-CK, F^l^K, F' = ^Ci, (2.22) 
m m m 

due to the fact that along a circular orbit r = const. The quantities a, F and F' 
are necessarily small. Although the quadrupolar terms / and /' are small only for 
a quasi-spherical body, the further rescaling by (k = {M/r^y^^ makes indeed them 
small in any case. In fact, the radius of the orbit is assumed to be large enough 
in comparison with certain natural length scales like \s\/m (also known as the M0ller 
radius [T5] of the body), (|/|/m)^/^, (|/'|/m)^/^ associated with the body itself in order 
to avoid backreaction effects. Furthermore, we also require that the characteristic 
length associated with the quadrupole are small if compared to the M0ller radius of 
the body, in order that the multipolar expansion of the body's stress-energy tensor 
underlying the Dixon's model be consistent. 
Eqs. ([2?T5l) and ([2T7l) then become 

= {vvp - vl)a + vk{v~ Vp) + ivji'MRF , (2.23) 

7 

and 

= {vvp - i^Di^F' + 2^lF) + il[2vvpvl + v^vi^p + vl) + 2vlvl]F 

- \llv{^''4 + H^^l - ^^^v^l - ^'^K + ^k) ■ (2.24) 

The above relations define the kinematical conditions allowing circular motion of the 
extended body taking into account its spinning and quadrupolar structures. The 
contribution due to the spinning structure disappears (i.e. remains arbitrary) when 
vvp — v^, i.e. if i^K is taken to be the geometrical mean of v and Vp. In fact, in this 
case, we have 

^ = -4^, (2.25) 

and 

0^{iy% + 2ul){v - Up) + {2^1 -y^- v],)vp . (2.26) 

After some manipulation using the condition i/Up = the latter equation turns out 
to be identically satisfied, whereas Eq. ()2.25|) becomes 



with F' arbitrary. 



Equations (|2.23p and (|2.24p can be examined in other special cases; for example: 
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= {jyi^p - vl){-F' + 2-ilF) + il[2uvyK + 4^^^p + ^l) + ^'^WpW 
9 

^ /22ir,22 o 2 22|4\ 

We notice that if = these equations imply v = Vp and 



= {vl - vl) 



that is Vp = -iiVK, or if ^ 



F' 



^ 2/2 



3 'PV P 

aUowing nongeodesic motion only due to the quadrupole moment tensor, 
(u) g 7^ 0, F ^ 0,FV . 

= {vi^p - v\)a + vk{v - Vp), 
2 



(2.28) 
(2.29) 

(2.30) 



= {yvp - Vi^)F' + -77p(i^ i^p + 2VpVj^ - ivVpVj^ - v Vj^ + v 



K) 



(2.31) 



(in) a ^0, F ^0,F' = . 

— {vup — vlli)<7 + vk{v — Vp) + iVK 



IpVp 



= jl[2vvp{l + 4) + v^pii^vp + 34) - 24]F 
9 

^ /22ir,22 o 2 22|4\ 

- o77p('^ + ^''p^K - ivvpvii -V Vk + ^k) ■ 



(2.32) 

The case of a spinning particle with vanishing quadrupole moment tensor, i.e. 
F = Q = F' has been already studied in 16J. 

Let us turn to Eqs. (|2.23p - (|2.24p and investigate the case of extended bodies 
with internal structure (dipolar and quadrupolar) compatible with a nearly geodesic 
motion. In the case of vanishing quadrupole moments Eqs. (12. 23^ - 1)2. 24|) reduce to 

m 

^ [vvp - v\)a + vk{v ~- Vp) , 



- ^ 2^.2^^ 



3z^^'„i^ 



■pi^K 



K ■ 



In the limit of small spin a we find 



Vp ^iy + 0{a^) 



(2.33) 



(2.34) 



If the contribution of quadrupolar terms can be considered negligible with respect 
to the dipolar ones and comparable with second order terms in the spin itself, one can 
consider corrections to Eq. (j2.34|) as given by 



3 3 
v± ^ ±vk ~ 2^^'^ ^ 8 ^ 



v± ± 3{F — a'^)vK , 



(2.35) 



where the signs ± correspond to co/counter rotating orbits. The corresponding 
angular velocity C,± — {Ck /vk)v± and its reciprocal are 



^-±^ + ^'^±^(5a^-2F) 

(,± (,K 2(,K OKK 



(2.36) 



Motion of extended bodies in Schwarzschild spacetime 



7 



Furthermore, the period of revolution around the central source will consist of three 
different terms 

T=^ = TK\l±Xd + \\. (2.37) 
K±l 

where 

Tk^^, Ad = ^a, X^^l{5a^-2F) . (2.38) 

A direct measurement of T will then allow to estimate the quantity F determining 
the quadrupolar structure of the body, if its spin is known. Note that the fraction 
due to the spin is different depending on whether the body is spinning up or down, 
whereas the term Xq due to the quadrupole has a definite sign once the shape of the 
body is known {F cannot change its sign). 

In the case of the Earth the nondimensional quantities (|2.22|) turn out to be given 
by cr « 2.3 X IQ-^^ and F = F' « -1.8 x 10"^°, since s/we « 3.4 x 10^ cm and 
f — f ~ — J2m©r^, with J2 w 10^^, and the distance between the Earth and the 
Sun is r « 1.5 X 10^'^ cm. The correction to the geodesic value Tk ~ 9.425 x 10^^ 
cm due to the spin is « 3.4 x lO"'^^, whereas the correction due to the quadrupole 
turns out to be A^ « 1.3 x 10"^'' as from Eq. (jOS)) . 

An interesting opportunity to test the quadrupole effect of extended body would 
arise from the motion of a binary pulsar system around Sgr A*, the supermassive 
{M ~ 10'^ Mq) black hole located at the Galactic Center [HI [Hj. To illustrate 
the order of magnitude of the effect, we may consider the binary pulsar system PSR 
J0737-3039 as orbiting Sgr A* at a distance of r ~ 10^ Km. The PSR J0737-3039 
system consists of two close neutron stars (their separation is only dAB ^ S x 10^ 
Km) of comparable masses uia ~ 1.4 Mq, itib — 1.2 Mq), but very different intrinsic 
spin period (23 ms of pulsar A vs 2.8 s of pulsar B) (T^]. Its orbital period is about 
2.4 hours, the smallest yet known for such an object. Since the intrinsic rotations 
are negligible with respect to the orbital period, we can treat the binary system as 
a single object with reduced mass (xab — 0.7 Mq and intrinsic rotation equal to the 
orbital period. The spin parameter thus turns out to be equal to cr « 6 x 10^^, 
whereas the quadrupolar parameters are F — F' k, 9.6 x 10~^°, since we have 
taken f = f = habc^ab ^ ^ rough estimate. The correction to the geodesic value 
Tk ~ 1.6 X 10^^ cm due to the spin is A^ ~ 9 x 10^*, whereas the correction due to 
the quadrupole turns out to be Ag « —7.2 x 10~^° as from Eq. (|2.38p . 

The black hole at the Galactic Center is actually expected to be rotating |20j . 
Therefore, a more detailed analysis would take into account the effect of rotation of 
Sgr A* on the estimate of the various contributions to the period of revolution of 
the PSR J0737-3039 binary system. We will discuss such an extension of the present 
treatment in a forthcoming paper. 



3. Quadrupolar corrections to the clock-efFect for spinning test particles 

Bini, de Felice and Geralico have investigated in |T6] the gravitomagnetic clock 
effect appearing for oppositely orbiting both spin-up or spin-down particles in the 
Schwarzschild spacetime. They found that spinning test particles move on circular 
orbits which, to first order in the spin parameter a, are close to a geodesic, with 

^ ±:^±-^k|, (3.1) 



C(±,±) (k Kk 
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where the signed magnitude cr = ±|cr| of the spin parameter has been introduced. The 
signs in front of 1/Ck correspond to co/counter-rotating orbits while the signs in front 
of \<7\ refer to a positive or negative spin direction along the z-axis; for instance, the 
quantity C(+,-) denotes the angular velocity of U corresponding to a co-rotating orbit 
(+) with spin-down (— ) alignment, etc. Therefore one can measure the difference in 
the arrival times after one complete revolution with respect to a static observer. The 
coordinate time difference in given by: 

At(+,+^_,+) = 27: (-^ + = , (3.2) 

VC(+,+) C(-,+)/ 

and analogously for At(+ _._ _). 

Let us analyze now the introduction of quadrupolar terms. Eq. (|2.36p implies 
that Eq. (|3.ip becomes 

-r^ = ± -^kl ± -^(5a2 - 2F) , (3.3) 

so that the last term does not contribute to the clock effect, since the ± signs in front 
of it correspond to co/counter-rotating orbits (like those in front of 1/(k), and thus 
cancels anyway. Therefore, no modifications arc induced by quadrupolar terms in Eq. 

(E21). 



4. Concluding remarks 

We have studied the motion of quadrupolar particles on a Schwarzschild background 
following Dixon's model. In the simplified situation of constant frame components 
(with respect to a natural orthonormal frame) of both the spin and the quadrupole 
tensor of the particle we have found the kinematical conditions to be imposed to the 
particle's structure in order the orbit of the particle itself be circular and confined on 
the equatorial plane. Co-rotating and counter-rotating particles result to have a non- 
symmetric speed in spite of the spherical symmetry of the background, due to their 
internal structure. This fact has been anticipated when studying spinning particles 
only, i.e. with vanishing quadrupole moments |16j . 

We have then discussed the modifications due to the quadrupole which could 
be eventually observed in experiments. Such experiments, however, cannot concern 
standard clock effects, because we have shown that there are no contributions arising 
from the quadrupolar structure of the body in this case. In contrast, the effect of 
the quadrupole terms could be important when considering the period of revolution 
of an extended body around the central source. Measuring the period will provide 
an estimate of the quantities F, F' determining the quadrupolar structure of the 
body, if its spin is known. On the other hand, the complete knowledge of the internal 
structure of the body will allow to estimate the period of revolution. In the latter 
case, the comparison between the measured period (known from observations) with 
the value predicted by the Dixon's model could provide a test for the model itself. 

It would be of great interest to extend this analysis to systems with 
varying quadrupolar structure and emitting gravitational waves without perturbing 
significantly the background spacetime. Usually variable quadrupole moment is 
generated in a test astronomical body of mass m because of tides produced by the 
central source of mass M ^ m. Now the net gravitational radiation due to motion of 
m is due to its orbit around M, the time varying tides and the interference between 
these two [2TJ[22]. We deserve such an investigation to future works. 
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